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Elliptic flow (1)2) and hexadecupole flow (174) of light clusters have been studied in details for 25 
MeV/nucleon 86 Kr + 124 Sn at large impact parameters by Quantum Molecular Dynamics model 
with different potential parameters. Four parameter sets which include soft or hard equation of 
state (EOS) with/without symmetry energy term are used. Both number-of-nucleon (A) scaling of 
the elliptic flow versus transverse momentum (p t ) and the scaling of V4/A 2 versus (pt/A) 2 have been 
demonstrated for the light clusters in all above calculation conditions. It was also found that the 
ratio of V4/V2 2 keeps a constant of 1/2 which is independent of p t for all the light fragments. By 
comparisons among different combinations of EOS and symmetry potential term, the results show 
that the above scaling behaviors are solid which do not depend the details of potential, while the 
strength of flows is sensitive to EOS and symmetry potential term. 
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Anisotropic flows are very useful to explore heavy- 
ion collision dynamics since it results from the transi- 
tion of the original space-time asymmetry into a mo- 
mentum space anisotrop y fo r the non-central collision 
[3, 1 i i, i, i 0, 1 1, MM, E 0, • Different mech- 
anisms will contribute the final momentum anisotropy, 
i.e. flow. Many studies of the dependence of the directed 
flow (t>i) and the elliptic flow (1*2) on beam energies, mass 
number, isospin and impact parameter have been carried 
out and much interesting physics has been demonstrated 
on the properties and origin of the collective motion in 
both nucleonic or partonic levels. Very recently, we car- 
ried out a Quantum Molecular Dynamics model calcu- 
lation with hard equation of state and symmetry energy 
term and found that there is a nucleon-number scaling for 
the elliptic flow of light particles up to the mass number 
A = 4 15]. In this work, we shall present more details for 
the nucleon number dependence of the anisotropic flows 
V2 and W4 for 86 Kr + 124 Sn collisions at 25 MeV/nucleon 
and large impact parameters (b = 7 — 10/m) with differ- 
ent EOS and symmetry energy interaction. The scaling 
behaviors look robust since they do not depend the pa- 
rameters used in the model, and the sensitivities of EOS 
and symmetry potential for the vi and W4 are discussed. 

Anisotropic flow is defined as the different n-th har- 
monic coefficient v n of the Fourier expansion for the par- 
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ticle invariant azimuthal distribution 

— — oc 1 + 2 v n cos(n0), 
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where <j) is the azimuthal angle between the transverse 
momentum of the particle and the reaction plane. Note 
that the z-axis is defined as the direction along the beam 
and the impact parameter axis is labelled as x-axis. 
The first harmonic coefficient V\ represents directed flow, 
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where p t -- y ^ , ry 

momentum. While the i>2 which measures the eccentric- 
ity of the particle distribution in the momentum space 
represents elliptic flow, 
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and i>4 represents the 4-th momentum anisotropy, namely 
hexadecupole flow: 
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The intermediate energy heavy-ion collision dynamics 
is complex since both mean field and nucleon-nucleon 
collisions are playing the competition roles. Furthermore, 
the isospin dependent role should be also incorporated for 
asymmetric reaction systems. Isospin dependent Quan- 
tum Molecular Dynamics model (IDQMD) has been affil- 
iated with isospin degrees of freedom with mean field and 
nucleon-nucleon collision [H, Q3, GJ, GJ, E3, EH HI H • 
The IDQMD model can explicitly represent the many 
body state of the system and principally contain corre- 
lation effects to all orders and all fluctuations, and can 
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describe the time evolution of the colliding system well. 
When the spatial distance (Ar) is closer than 3.5 fm 
and the momentum difference (Ap) is smaller than 300 
MeV/c between two nucleons, two nucleons can coalesce 
into a cluster [l6j]. With this simple coalescence mech- 
anism which has been extensively applied in transport 
theory, different size clusters can be recognized. 

In the model the nuclear mean-field potential is pa- 
rameterized as 

U(p, t z ) = a(^) + (3(-^y + hi r z )V c 
Po Po 2 
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where po is the normal nuclear matter density 
(0.16/m~ 3 ), p n , p p and p are the neutron, proton and 
total densities, respectively. r z is z-th component of the 
isospin degree of freedom, which equals 1 or -1 for neu- 
trons or protons, respectively. The coefficients a, f3 and 7 
are parameters for nuclear equation of state. C sym is the 
symmetry energy strength due to the density difference 
of neutrons and protons in nuclear medium, which is im- 
portant for asymmetry nuclear matter [24], [25], [26| (here 



32 MeV is used to consider symmetry energy 



effect or isospin-dependent potential, and C sy m = for 
no symmetry energy effect or isospin-independent poten- 
tial). V c is the Coulomb potential and JJ Yuk is Yukawa 
(surface) potential. In this work, we take a = 124 MeV, (3 
= 70.5 MeV and 7 = 2 which corresponds to the so-called 
hard EOS with an incompressibility of K = 380 MeV, 
and a = -356 MeV, (3 = 303 MeV and 7 = 7/6 which 
corresponds to the so-called soft EOS with an incom- 
pressibility of K = 200 MeV. In the present study, four 
combinations with different potential parameters, i.e. pa- 
rameters of hard or soft EOS with or without symmetry 
energy effect (i.e. C sym = 32 or MeV), for the collision 
system of 86 Kr + 124 Sn at 25 MeV/nucleon with impact 
parameter from 7 fm to 10 fm were carried out. The 
physics results were extracted at the time of 200 fm/c 
when the system has been in the freeze-out stage. 

The Fig.l (a), (b), (e) and (f) shows transverse momen- 
tum dependence of elliptic flows for mid-rapidity light 
fragments in four different calculation conditions: (a) for 
soft EOS with symmetry potential (softJso); (b) for 
hard EOS with symmetry potential (hardjiso); (e) for 
soft EOS without symmetry potential (softjniso) and (f) 
for hard EOS without symmetry potential (hardjniso). 
In all cases, elliptic flow is positive and it increases with 
the increasing p t , which is apparently similar to RHIC's 
results [HI, El]. Of course, the mechanism is very differ- 
ent. In intermediate energy domain, collective rotation 
is one of the main mechanisms to induce the positive el- 
liptic flow 0, H, H3, [H, HI] . However, at RHIC energies 
it is the strong pressure which is built in early initial 
almond anisotropy of the geometrical overlap zone be- 
tween both colliding nuclei that drives the positive ellip- 
tic flow fl2l ]. The corresponding nucleon- number scaled 
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FIG. 1: (a), (b), (e) and (f): Elliptic flow as a function of trans- 
verse momentum (pt) for the simulation with different parameters 
of EOS with or without symmetry energy term, (a) for soft EOS 
with symmetry potential (softjso); (b) for hard EOS with sym- 
metry potential (hardSso); (e) for soft EOS without symmetry 
potential (softjniso); and (f) for hard EOS without symmetry po- 
tential (hardjniso). Squares represent for neutrons, circles for pro- 
tons, triangles for fragments of A=2, diamonds for A=3 and stars 
for A=4. Fig.(c), (d), (g) and (h) presents nucleon- number normal- 
ized elliptic flow as a function of transverse momentum per nucleon 
corresponding to the case of (a), (b), (e) and (f), respectively. 



elliptic flows are plotted in Fig.l (c), (d), (g) and (h) as 
a function of transverse momentum per nucleon. From 
these panels, it seems that the number of nucleon scal- 
ing for elliptic flow exists for light fragments at low Pt/A 
(pt/A < 0.2GeV/c). This behavior is apparently simi- 
lar to the number of constituent quarks scaling of elliptic 
flow versus transverse momentum per constituent quark 
(j>t/n) for different mesons and baryons which was ob- 
served at RHIC [l2j . Since all calculations show the sim- 
ilar scaling behavior, this scaling behavior is robust, and 
it is independent of the details of EOS and symmetry 
potential. 

To quantitatively look the difference of the flows in 
different calculation conditions, we compare the values 
of V2/A for the four simulation conditions (see Fig. 2). 
The figures show that the difference between different 
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FIG. 2: Comparisons of the values of V2 /A versus pt/A in the 
different simulation conditions. The meanings of the symbols are 
depicted in right bottom corner. 



simulations is big for neutrons and protons but a little 
small for the fragments of A = 2, A = 3 and A = 4. 
The reason is that the emitted protons and neutrons can 
feel the role of mean field (EOS) directly, while the light 
fragments have weak sensitivity since they are indirected 
products by the coalescence mechanism in the present 
model. Approximately at the samep t /A, the elliptic flow 
is larger for soft EOS than the one for hard EOS, and 
it is larger for EOS with symmetry potential than the 
case without symmetry potential. Considering that the 
symmetry potential is basically positive for the studied 
reaction system (more neutrons than protons) , symmetry 
potential will make the whole EOS stiffer. In this case, 
we can say that the stiffer the EOS, the smaller the flow. 
In other words, we can say the strength of elliptic flow per 
nucleon is sensitive to the EOS and symmetry potential. 

So far, there is rare studies about higher order flows, 
such as U4, experimentally and theoretically in this en- 
ergy domain. Here we try to explore the behavior of v^. 
First we draw v^/A as a function of Pt/A mimicing the 
behavior of elliptic flow (see (a), (b), (e) and (f) of Fig. 3) 
for four different calculation conditions. It shows that 
Vi/Ais positive and increases with p t /A, but there seems 
no simple scaling behavior as v-i shows. Considering that 
RHIC experimental data have demonstrated that a scal- 
ing relation among hadron anisotropic flows holds, i.e., 

v n {pt) ~ v 2^ 2 (Pt) [33] > we pl°t v a/A 2 as a function of 
(pt/A) 2 in Fig. 3 (c), (d), (g) and (h) for the correspond- 
ing calculation conditions of Fig. 3(a), (b), (e) and (f). 
Now the points of different light fragments nearly merge 
together at low (p t /A) 2 , which means a certain of scal- 
ing law holds between two variables. All the calculation 
cases show that there is the scaling behavior for V4/A 2 
versus (pt/A) 2 , and this behavior is robust regardless the 
parameters which we used for EOS. 

Since the above scaling behavior assumes v n (p t ) ~ 

i>2^ 2 (pt)' so we pl°t v 4/ v 2 as a function of p t in Fig. 4 
for the four simulations. The figures show that the ratios 
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FIG. 3: Same as Fig.l but for V4/A versus pt [(a), (b), (e) and (f)] 
and Vi/A 2 versus (pt/A) 2 [ (c), (d), (g) and (h)]. 
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FIG. 4: The ratios of vn/v 2 for neutrons (squares), proton (circles), 
the fragments of A = 2 (triangles), A = 3 (diamonds) and A = 4 
(stars) versus pt for the four simulations in different calculation 
conditions. 
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FIG. 5: Comparison of the values of V4/A versus pt /A. The mean- 
ings of the symbols are depicted in right bottom corner. 



of V4/V 2 for different fragments up to A ~ 4 are about a 
constant of 1/2 in all simulation cases. Because V2/A can 
be scaled with pt/A, V4/A 2 should scale versus (pt/A) 2 , 
which is exactly what we see in Fig. 4. One point is 
worth to be mentioned comparing to the RHIC studies 
where the data shows V4/V 2 ~ 1.2 :3(|, vn/v 2 ~ 1/2 
for the light nuclear fragments in this nuclconic level co- 
alescence mechanism rather than the value of 3/4 for 
mesons or 2/3 for baryons in quark coalescence model 
[3lj . Coincidentally, the predicted value of the ratio of 
V4/V2 for hadrons is also 1/2 if the matter produced in 
ultra-relativistic heavy ion collisions reaches to thermal 
equilibrium and its subsequent evolution follows the laws 
of ideal fluid dynamics 32]. It is interesting to note the 
same ratio was predicted in two different models at very 
different energies, which is of course worth to be further 
investigated in near future. One possible interpretation 



is that the big nucleon-nucleon cross sections in low en- 
ergy HIC make the system to reach thermal equilibrium 
and may induce the fluid-like behavior of nuclear medium 
before the light fragments are coalesced by nucleons. In 
this case, the value of v^/v 1 of light fragments could be 
- 1/2 as Ref. [13] shows. 

The values of V4/A versus pt/A with different simula- 
tion parameters are also presented for light fragments, see 
Fig. 5. The figures are similar to those in Fig. 2, and the 
effects of EOS and symmetry potential on V4/A are also 
similar to their effects on V2- However, comparing with 
the V2 J s sensitivity to the EOS and symmetry potential, 
V4/A is not so salient. 

To summarize, we investigated the behavior of 
anisotropic flows as a function of transverse momentum 
for light fragments for the simulations of 25 MeV/nucleon 
8G Kr+ 124 Sn collisions in peripheral collisions by IDQMD 
model in the potential parameters of hard or soft EOS 
with or without symmetry energy term. It was found 
that for all the four type simulations V2 and V4 of light 
fragments are positive and increase with pt/A. When 
we plot V2 per nucleon fa/ A) versus Pt/A for all light 
particles, all curves collapse onto the same curve. Sim- 
ilarly, the values of V4/A 2 merge together as a function 
of (pt/A) 2 for all light particles. Furthermore, it was 
found that V4 can be well scaled by and the value of 
V4/V 2 ~ 1/2 which does not depend on transverse mo- 
mentum. The above scaling behaviors can be seen as an 
outcome of the nucleonic coalescence, and it illustrates 
that the number-of-nucleon scaling for elliptic flow exists 
in intermediate energy heavy ion collision. In addition, 
the values of V2/A and V4/A were compared in different 
simulation conditions, and it was shown that the values 
of the V2 are sensitive to the EOS and symmetry poten- 
tial, especially for neutrons and protons. 
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